Physical region for three-neutrino mixing angles 
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We derive a set of symmetry relations for the three-neutrino mixing angles, including the MSW 
matter effect. Though interesting in their own right, these relations are used to choose the physical 
region of the mixing angles such that oscillations are parameterized completely and uniquely. We 
propose that the preferred way of setting the bounds on the mixing angles should be #12 £ [0, tt/2], 
613 £ [— 7r/2, tt/2], #23 £ [0, 7r/2], and 8 £ [0,7r). No CP violation then results simply from setting 
5 = 0. In the presence of the MSW effect, this choice of bounds is a new result. Since the size of 
the asymmetry about #13 = is dependent on the details of the data analysis and is a part of the 
results of the analysis, we argue that the negative values of #13 should not be ignored. 
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An abundance of data now demonstrates that neutri- 
nos oscillate between flavor states. Analyses of the data 
in the context of three-neutrino mixing can be found in 

H7I Ha | . Conventionally, the standard model of the elec- 
troweak interaction is extended to include a mass term for 
the neutrinos and a unitary mixing matrix which relates 
the flavor states of the neutrinos to the mass eigenstates. 
We here utilize an algebraic formalism to derive sym- 
metries of the mixing matrix which leave the predicted 
oscillation probabilities invariant. In particular, we are 
interested in understanding the bounds on the mixing 
angles which parameterize the mixing matrix. We derive 
the known bounds for the general case which includes 
CP violation. For the case of no CP violation, we pro- 
pose a new result, for the case in which the MSW ef- 
fect is present, of choosing the bounds as 612 E [0,7r/2], 
013 € [-tt/2, tt/2], 6> 23 e [0,7r/2], (5 = 0. We argue that 
this choice has distinct advantages for the phenomcno- 
logical analysis of neutrino oscillation data. 

The transformation from the mass states to the fla- 
vor states is represented by a unitary matrix U a j, v a = 
U a ji>j, where we use Roman letters, j, for mass eigen- 
states and Greek letters, a, for flavor states. For two neu- 
trinos, this mixing matrix reduces to an element of the 
commutative group U(l). For vacuum oscillations, the 
mixing angles lie between and 7r/4 in order to describe 
all discernible scenarios. For propagation in matter, this 
range on the mixing angle must be expanded to [0, 7r/2]. 

In a three-neutrino theory, the mixing matrix is an el- 
ement of a non-commutative group. Symmetries in the 
parameterization are then less obvious, especially for neu- 
trinos propagating in matter. Transcendental expressions 
for locally defined mass eigenstates exist , but the un- 
derlying symmetries of these states are not transparent. 
As such, we extend the algebraic formulation of vacuum 
neutrino oscillations developed in [^j to include matter 
effects. Thus we may derive symmetries among the mix- 
ing angles and the bounds on the mixing angles which 



these symmetries imply. The bounds on the mixing an- 
gles are well known in vacuum; for this case, the neutrino 
oscillation problem is then mathematically equivalent to 
the quark mixing problem. 

The limits on mixing angles, including the propagation 
through matter, have previously been considered in |2l| . 
Our results are in agreement with theirs; however, our 
use of the symmetry relations permits us to arrive at 
additional conclusions. A discussion of the four neutrino 
case can be found in Ref. |22^ . 

In [2(J , the vacuum oscillation probability, valid for an 
arbitrary fixed number of neutrino mass eigenstates, is 
given by 



■Ht pa g—iHt 



(i) 



where P a is the flavor projection operator given by 
(P a )jk — (U^Q a U)jk — U*jU a k in the mass-eigenstate 
basis, P + projects onto the positive-energy states, and 
H is an rt-particle Dirac Hamiltonian for masses rrij . We 
note that adding multiples of the identity to the Hamilto- 
nian leaves Eq. Q invariant, and we recall that the trace 
is cyclical. For the case of three neutrinos, the relativistic 
limit, p 3> m, of the above equation is equivalent to the 
usual oscillation formula 



V a ^p(L/E) = ti[Ue lML ' 2E U^Q a Ue- lML / 2E U^Q 13 ] 



(2) 
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(3) 



where A4 is the diagonal matrix with entries 
(m 2 ,rri2,m 2 ) and ipjk '■= AjkL/4E with Ajk := m 2 —m 2 . 
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The generalization to neutrinos propagating in mat- 
ter requires that the generator of time (space) transla- 
tions be modified to include a flavor-dependent potential 
which accounts for interactions with the electrons in mat- 
ter [2^]. The neutral current interaction does not effect 
oscillations and is neglected. The charged current inter- 
action involves only the electron flavor. As such, we must 
add to the Hamiltonian an effective potential which op- 
erates exclusively on the electron flavor. Explicitly this 
potential is A{x) = V2 G E p{x)/m n , with p{x) the elec- 
tron density at position x, G the weak coupling constant, 
m n the nucleon mass, and E the neutrino energy. In the 
mass basis, the dynamical equation in the presence of 
matter then becomes 



idtv m = (H + V)v m = Hl/r, 



(4) 



where V is given by A(x)P e with P e the electron flavor 
projection. Given nontrivial mixing among the mass and 
flavor states, it is clear that the free Hamiltonian H does 
not commute with the potential V. 

In the relativistic limit, one has x ~ t so that temporal 
derivatives are, in effect, spatial derivatives. Formally, 
the solution to the differential equation in 10} is 



v m (Y) = exp 



Hdx \ f m (0), 



(5) 



where T parameterizes the oriented path of the neutrino 
through the matter with starting point x = 0. As the 
Hamiltonian is hermitian, the adjoint of the above oper- 
ator is simply exp (i J r Hdx^ which merely reverses the 
orientation of the path T. 

Written in a manner akin to that in J2J , the expression 
for the oscillation probability in matter is 



P a -H3(F,E) = tr 



exp < i 



exp 



(6) 



In the relativistic limit, the argument of the path ordered 
exponentials in this trace can be somewhat simplified 

J UHUUx = UMU^L/2E + J A(x)dxQ e , (7) 



modulo a multiple of the identity; as noted above, mul- 
tiples of the identity have no physical consequence with 
respect to oscillations. 

We use the standard parametrization of the three- 
neutrino mixing matrix |24| 



C12C13 S12C13 S136 lS 



U (623, 6x3,6x2, 5) = I -S12C23 - CX2S2SSXS& C12C23 - Sl2S23Sl3e S23C13 I , (8) 
S12S23 - Ci2C 2 3Sl 3 e l<5 -C12S23 - Sl2C23Sl3e l<5 C23C13 



where Cjk — cosOjk, Sjk = shxOjk, and Qjf., 5 are 
real. We deem equivalent two quadruples of parameters 
(#23, $137 6\2, S) if the oscillation probabilities P a ^,p(T, E) 
are identical for all values of a and /3 over some neutrino 
path r. This equivalence will be expressed via the rela- 
tion =. Central to our method is the fact that the trace 
of an operator is invariant under conjugation by any uni- 
tary U, 

tx[UAU^] = it[A] . (9) 

As the flavor projections Q a,/3 in JJ5J must remain un- 
changed, they must commute with the conjugating ma- 
trix. 

The group SO(3) has three generators. It is conve- 
nient to express the mixing matrix in terms of the expo- 
nentiated generators Uj(6), a rotation by angle 6 about 
the jth axis in K 3 . To account for the CP phase, we 
make the additional definition for Sg £ SU(3) by setting 
S 5 = diag(e- i5 / 2 , 1, e lS/2 ) and letting 

Ui(6) = S s U 2 (9)Sl . (10) 



The neutrino mixing matrix JSJ) may be written as 

U(d23, 6x3, 6x2, 5) = Ux(6 2 3)Ui(6x3)U3(6x 2 ) ■ (If) 

From Eq. (|I0|) . we may extract our first equivalence 
relation. For the CP phase S — it, one finds 

SM6)Sl = U 2 [-6) ■ (12) 

hence, one has 

U s +*{6) = U b 2 {-6) . (13) 

The addition of tt to the CP phase can accommodate a 
change in sign of 6x3 

{623,6x3, 6x2,6) = {623, -6x3, 6x2, S + tt) . (14) 

In order to proceed to further symmetries, we employ 

U^U^U^ =U k (-6) (15) 

for j =/= k; this is easily confirmed by explicit calculation. 
Additionally, as Uj(w) is diagonal, it commutes with S$ 
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so that 

UjWUiWUj^ = Ul(-9) (16) 

for j = 1,3. 

For the same reason, Uj(w) commutes with the mass- 
squared matrix M and the projections Q a . The remain- 
ing step is to determine the action of Uj (tt) on the gen- 
erator of spatial translations. Referring to Eq. J7J), one 
sees that Uj(w) commutes with the MSW potential and 
thus only affects the term containing the mass-squared 
matrix M.. Combining these facts with Eq. ©, we see 
that the mixing matrices U and Uj(iryUUj(iT) yield iden- 
tical oscillation probabilities. Using (|15|1 and 1)16(1 . we 
may express this equivalence in terms of the parameters 
themselves 

(023,013,012,*) = (-023,-013,012,*) (17a) 
= (-023,013,-012,*) (17b) 
= (023,-013,-012,*) ■ (17C) 

From the same equations, we may determine the effect 
of adding tt to a mixing angle. We find the following 
equivalence amongst parameters 

(023,013,012,*) = (023 +7T, 013,012,*) (18a) 
= (-023,013 + 77,^12,*) (18b) 
= (023, 013, 012 +7T, S) . (18C) 

Trivially, the mixing angles satisfy a 2tt periodicity; 
however, from the relations in H18aH18c)) . it is clear that, 
without loss of generality, one may further restrict all 
djk to lie |25| in the interval [0, tt). In the case of #13, we 
are able to further narrow these bounds to [0,7r/2]. An 
application of relation l(17a(l followed by l(18b|l yields 

(023, 013, 012, *) = (023, - 013, 012, *) . (19) 

In general, given a real number x, the map x <— > a — x is 
a reflection about the point a/2 on the real line. Hence, 
we see that if #13 should lie between it/2 and tt, then 
the above relation shows that we have an equivalent os- 
cillation probability for an angle reflected about tt/2. In 
short, we may choose #13 to lie in the first quadrant. 

For the remaining mixing angles, we may apply rela- 
tions (jlTbf) - (|18a|l . and IJlScjl to achieve 

(023, 013, 6 12 ,S) = (TT - 023, 013, TT 012, *) • (20) 

From this, one notes that a reflection about tt/2 in both 
023 and 0i2 results in an equivalent oscillation probability. 
These reflections cannot be performed independently and 
still result in an equivalent theory. Hence, we only have 
the freedom to restrict either 823 or 6*12 to the interval 
[0,7r/2], but not both. 

Relation (|20|) is valid for fixed 5. Should we relax this 
condition on <5, then we have the liberty to restrict both 



mixing angles to lie between and 7r/2. Using iQ , (|17afl . 
and l|18a(l . we have 

(023, 013, 012, *) = (7T - 23 , 013, 012, * + 7r) . (21) 

From similar logic, one can also deduce 

(023, 013, 012, *) = (023, 013, TT - 012, * + 7l) . (22) 

Permitting a change in the CP phase results in the ability 
to independently reflect these two mixing ang les about 
7r/2; hence, we are guaranteed, as noted in |2lj,|26j, that 
all 0jk can be restricted to the interval [0, tt/2] if one al- 
lows the full range on the CP phase S S [0, 27r). This 
is the common understanding for the CKM mixing ma- 
trix for quarks and for neutrino oscillations including CP 
violation. 

We propose an alternate but equivalent set of bounds 
for the CP violating case. As a consequence of relation 
|Q , we see that the full range of the CP phase may 
be replaced by S € [0, tt) by allowing a negative mixing 
angle 0i3 G [— 7r/2, 7r/2] with the two other mixing angles 
remaining in the first quadrant. An alternative would 
be to utilize Eq. (j2U|l and allow either i2 £ [0,7r) or 
023 G [0, tt) with the CP phase in the first and second 
quadrants; however, given the smallness of 0i3 in most 
models we find the extension to negative values of 0i3 
most natural. 

Given that there is no indication of CP violation in 
neutrino oscillations, we give this case special consider- 
ation. The conventional view of taking the the mixing 
angles 0jk € [0, 7r/2) requires that both 6 = and 6 = tt 
must be included to cover the total allowed space for no 
CP violation, as was shown in [ilj . A common oversight 
has been to drop the 5 = tt branch. For no CP violation it 
is more natural to have 0i2 6 [0, 7r/2], i3 G [— tt/2, 7r/2], 
023 £ [0, 7r/2], with only the <5 = branch. This choice 
is motivated by the smallness of 0i3. We note that pro- 
viding an upper limit on the value of sin 2 0i3, a common 
practice |2^, is insufficient since 0i3 < is physical and 
the mixing is not symmetric in 0i3. A two neutrino anal- 
ysis cannot distinguish the sign of 0i3, but a linear term 
in 0i3 is present for three neutrinos and, depending on 
the details of the analysis, may be significant. 

Experiments are not yet sensitive enough to determine 
the Dirac phase 5. Until there exists data to indicate a 
value of 5 other than or tt in the convention of Ref. [2l| , 
or other than in the convention proposed here, the 
special case of no CP violation is the default assump- 
tion. An example of the necessity, in principle, of incor- 
porating the negative 0i3 region arises from examining 
the question of mass hierarchy. To demonstrate this, we 
adopt the convention for ordering the mass eigenstates 
by increasing mass, ni\ < rri2 < m^. This, together 
with the bounds chosen for the mixing angles, defines 
a unique basis for an analysis, i.e. each possible phys- 
ical solution is included once and only once. Let us 
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take as a typical set of parameters derived from oscil- 
lation data, A 2 i = 7 x l(T 5 eV 2 , A 32 = 3 x lCT 3 eV 2 , 
012 = 0.60, 0i3 = 0.10, and 23 = 0.80. This set 
of mass-squared differences is an example of the regu- 
lar hierarchy. It is known that there is set of param- 
eters with an inverted mass hierarchy that yield oscil- 
lation probabilities which are nearly equivalent. Using 
Eq. (28) in j^], we may calculate these nearly equiva- 
lent parameters A 2i = 3 x 10~ 3 eV 2 , A 32 = 7 x 10~ 5 eV 2 , 

012 = 1.45, 0i3 = -0.60, and 23 = 0.70. We note the 
value of 0i 3 is negative and of relatively large magnitude. 
This is a consequence of choosing a fixed (ascending) 
mass order. Though we concede the common conven- 
tion for the inverted hierarchy is to reorder the masses, 
the point being made here is that the resolution of the 
mass hierarchy question is independent of arbitrary mass 
ordering schemes. Utilizing a single basis, as is usually 
done in quantum mechanics, requires negative values of 

01 3 (or the separate branch that arises from 6 — it) to be 
mathematically complete and correct. 

Vacuum oscillations are a special case of oscillations 
in matter. Consequently, they automatically satisfy the 
above relations. The remaining question is to determine 
whether there exist any additional symmetries for the 
vacuum case, as is in the two-neutrino theory. For three 
neutrinos, there are none. In our above analysis, we note 
that the operator Uj (it) commutes with the MSW poten- 
tial in Eq. J7J. As such, the presence of this potential is 
inconsequential in the derivation of the symmetries that 
follow; hence, the special case of vanishing electron den- 
sity admits no additional symmetries. To ensure that 
we have not overlooked a symmetry which could fur- 
ther reduce the bounds on the mixing angles, we have 
checked numerically that the mixing probabilities are in- 
deed unique over the regions given. 

In summary, we have put forth a new method for de- 
riving symmetries of the three-neutrino mixing matrix. 
These symmetries are of interest themselves. They also 
provide an elegant way to determine the physical bounds 
on the mixing angles in the presence of the MSW effect. 
In particular, we note that a convenient way of setting 
these bounds is to take 0i 3 € [—ir/2, ir/2], S 6 [0, it), and 
the other two angles in the first quadrant. In the case 
of no CP violation, where some amount of confusion ap- 
pears in the literature, this reduces to setting 6 = 0, and 
taking 0i 3 £ [— tt/2, tt/2] with the other two angles in the 
first quadrant. We note that in an expansion in terms of 
the ratio of the small mass square difference to the large 
mass square difference, the lowest term is dependent only 
on sin 2 0i 3 . In this limit, negative 0i 3 would be redun- 
dant. However, the leading correction is linear in 0i 3 . 
We examine this question in Ref. and find that, for 
reasonable parameters, the x 2 space may have a notice- 



able asymmetry about 0i 3 — 0, thus further motivating 
that the parameter space for the mixing angles be taken 
as derived here, a result that in the presence of the MSW 
effect, has not previously appeared in the literature. 
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